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Abstract 


A computer program is developed to predict the performance of a thermoacoustic 
generator. The defining relationships are based on equations derived using short stack 
and boundary layer approximations. The engine modeled is the five inch engine 
currently in operaticn at Los Alamos National Laboratory. The working fluid in the 
model is a helium-argon gas mtx, chosen to provide flexibility in placing the 
trausuucer relative to the thermoacoustic stack assembly while maintaining 
resonance near the desired operating frequency of 120 Hz. The transducer model is 
based on a linear alternator design proposed for use with the L. A. N. L. engine in the 
future. The thermodynamic properties of the gas mix are determined using virial 
equations and binary gas mixture relationships. The defining relations are solved 
simultaneously and an iterative process is used to optimize the gas mix. The variables 
for which the program then solves are pressure amplitude. hot heat exchanger 
temperaiure, and angular frequency as functions of input heat flux, system geometry, 
and load characteristics. From these, the power output and total system efficiency are 


determined. 
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I. INTRODUCTION 


A. OVERVIEW 

A great deal of work has been done in the development of practical 
thermoacoustic refrigerators which use gas as a working fluid. Examples include the 
Thermoacoustic Refrigerator (Hofler, 1986) and the Space ThermoacousStic Refrigerator 
(Susalla, 1988). Work on thermoacoustic prime movers which use a gas as a working 
fluid consists primarily of enclosing the thermoacoustic engine in a tube with two rigid 
ends and measuring the performance of the engine under various pressure, temperature 
and geometry conditions. Relatively little work has been done In developing a 
thermoacoustic generator. 

The history of thermoacuustics is long, but significant advances in practicai 
thermoacoustic devices are 1elatively recent. As early as 1777, experiments were 
performed by Byron Higgins, in which acoustic oscillations were excited in a large pipe 
by suitable placement oi an internal hydrogen flame (Putnam and Dennis, 1956). 
Higgens’ work eventually led to the modern science of pulse combustion, whose 
applications have included the German V-1 rocket used in World War II and the 
residential pulse-combustion turnace intruduced by Lennox, Inc. in 1982 (Swift.1982). 

The earliest antecedent of the thermoacousStic primrne movers in use today is the 
Sondhauss tube. Over 100 years ago, it was observed that a hot glass bulb attached toa 
cool glass tube sometimes emitted sound at the tip of the tube. Sondhauss 
quantitatively investigated the relationship between the pitch of the sound emitted 
and the dimensions of the tube. Lord Rayleigh explained the Sondhauss tube 
qualitatively tn 1896, but a quantitative theoretical description of these phenomena 


was not achieved until about 20 years ago when Rott and others published a series of 


papers (Swill, 1988). Rott's theory combines basic principles from physics, 





thermodynamics, and acoustics to describe the effects found in both thermoacoustic 
prime movers and heat pumps. 

One of the most significant advances in experimental thermoacouslics was the 
realization, by Carter, White, and Steele of North American Aviation, Inc. in 1962, that 
placing structures (such as a stack of plates) inside a Sondhauss tube could improve the 
performance of the tube in generating sound (Swift, 1988). 

In recent years, a number of people including Wheatley, Hofler, Swift, and 
Migliori, performed a series of experiments at Los Alamos National Laboratory 
(L.A.N.L.) in New Mexico. The purpose of the experiments was to investigate the 
thermoacoustic effects in both prime movers and heat pumps, and compare these 
experimental results to Rott's theory (Susalla, 1988). 

The first acoustic heat pump built at L.A.N.L. consisted of a speaker mounted at 
one end of a closed tube which contained a stack of fiberglass plates positioned near the 
closed end. They found that a temperature difference of 100 degrees centigrade across 
the stack of plates was easily produced in about a minute with the speaker driving at 
acoustic resonance (Wheatly, Swift. and Mighori, 1986). 

Shortly after the first heat pump was built at L.A.N.L,. Hofler built a device, 
which has come to be known as the Hofler tube, to use as a demonstration for his Ph.D. 
candidacy commuttee at the University of California, San Diego. The device (Fig. 1) 
preduces an impressive level of sound from the cren end, about 100dB at one meter 
(Swift, 1988). 

Other devices built and studied at L.A.N.L. include Hofler's refrigerator, the heat- 
driven acoustic cooler (Fig. 2), a quid sodium acoustic engine, and a liquid propylene 
engine (Wheatly. Swift. and Migliori, 1986). Currently being studied at L.A.N.L. 1s a five 


inch thermoacoustic prime mover, which will be described in detail later. It is the 


prime mover upon which the calculations described herin are based. 
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Figure 1: The Hofler tube. Loud sound is produced when a sufficient temperature 
gradient is established across the stack (Switt, 1988). 


B. MOTIVATION 


The results of these early experiments on thermoacoustic refrigerators at L-A.N.L, 
indicated a system performance that was far less than expecied (Holler, 1986). These 
results led Hofler to apply Rott’s theory to the expemmental systems, solving the 
resulting equations numencally. He also built the thermoacoustic refrigerator, made 
accurate measurements of its thermodynamic efficiency. and commpaied these resulis 
to those predicted by theory. The work which followed at the Naval Postgraduate 
School by Susala and others was done to unprove the overall efficiency of the system 
and make it suitable tor space applications (Susalla. 1986.. 


Swit s research group at L_A.N.L. is currently persuing a simular Line of 
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Figure 2: The heat driven acousuc cooler. High-temperature heat causes low 
temperature refrigeration (Swit, 1988) 





investigation for a thermoacoustic generator. A five inch diameter, four meter long 
prime mover has been built and operated in a rigid-rigid resonator. The system is 
operated at a mean pressure of twenty atmospheres with helium as the working fluid. 
Acoustic pressures on the order of ten percent of the mean pressure have been observed. 
Acoustic ower generated ts on the order of one kilowatt. 

The next phase of the develo, aent of an acoustic generator, using the five inch 
engine, is to obtain a transducer able to convert the acoustic energy to electrical energy. 
Two contractors, Mechanical Technology Inc. of Latham, New York and Sun Power Inc. 
of Athens, Ohio have been requested to provide quotations for delivery of a linear 
alternator for this purpose. The cost of such a device has been estimated at between 
$25,000 (Sunpower) and $115,000 (MTI). Because of the high cost (and current lack of 
funding) it was decided thal a computer program should be developed in the interim, 
which would predict the performance of the combined prime mover-transducer system. 
In addition, such a program, once experimentally validated, could be used to determine 
the effect of changing the geometry or materia: of the stack assembly, the position of 
the stack assembly in the resonator, or the characteristics of the mechanical or 
electrical load on the performance of the system as a whole. Such a program could 


result in. a great deal of savings in material and man hours when design modifications 


are contemplated. 





IL. THEORY 


A THERMODYNAMICS 

In designing any heat engine, compromises are made between thermodynamic 
ideals, and engineering and marketing realities. In designing any device. a practical 
balance is sought between the desire for high effictency, power, and reliability and the 
necessity of achieving low cost, small size and maximum simplicity. Thermoacoustic 
engines are extremely simple devices, as will be seen, and promise to deliver a 
substantial fraction of the ideal thermoacoustic efficiency. 


In any basic thermodynamics course we learn that there are two principle types of 


heat engines: heat pumps and prime movers. A heat pump removes heat Qc from a low 
termperalure reservoir al temperature Tc, receives work W from the surroundings, and 
delivers heat QH to a high ternperature sink at temperature TH. 

The process is reversed in a prime mover. Heat QH ts received from a high 
temperature reservoir at temperature TH, work W is done on the surroundings, and the 


remaining heat Qc is dumped to a low temperature sink at temperature Tc. 


The first and second laws of thermodynamics place an upper ltmit on the 
efficiency of any heat engine. This limit is best understood by applying the Carnot cycle 
to the process. The Camot cycle was first introduced in 1824 by N. L. Sad! Carnot (Sears 
and Salinger, 1975). The Carnot cycle (Fig. 4) consists of two reversible isothermal 
processes and two reversible adiabatic processes. The direction of the cycle (i.e. 
clockwise or counter-clockwise) determines whether the engine is a prime mover ora 
heat pump. Since the work in this thesis deals solely with prime movers, the remainder 


of the discussion tn this section will focus on them. 





From the first law of thermodynamics we know that energy must be conserved in 


the system. It follows that 
QH - Qc =W. (1} 
The second jaw requires that the total entropy in the system must increase or remain 


constant. The change in entropy for isothermal heat transport to/from a resorvoir is 


given by 


Prime mover Heat pump 


Figure 3: Thermodynamic cycles of a heat pump and prime mover. 








AS = G/T. (2) 


For a prime mover the second law may be written as 
ASc - ASH = Qc/TC - QH/TH 2 0. (3) 


The efficiency of interest tn a prime mover (the work output divided by the heat input) is 


N= W/QH. Combining Eqs. 1 and 3 it can be shown that 


T| = W/QH $ (TH - Tc)/TH (4) 


The temperature ratio on the right side of the above equation is known as Carnot’s 


efficiency ((c = (TH - Tc)/TH) and is the highest efficiency that a prime mover can 


achieve. 


The Camot cycle and Salinger. 











B $$THERMOACOUSTIC THEORY 
As stated in Chapter J, thermoacoustic theory was developed in detail by Rott and 
others in the late 1960's and 1970's. Since then, a great deal of work has been done by 
Swift, Hofler, Wheatley, and others in applying Rott's theory to experimental devices. 
The work contained in this thesis relies heavily on the relationships developed by 
applying the short stack and boundary layer approximations to the most general of 
thermoacoustic derivations (Swift, 1988). 
L A Qualitative Picture 
Consider a gas parcel in thermal contact with and oscillating along a plate at 
the acoustic frequency (Fig. 5). As the gas parcel moves back and forth, it will experience 
changes in pressure, volume. and temperature, The temperature change is due to both 
the adiabatic compression and expansion of the gas by the acoustic pressure and by the 
local temperature of the plate. 
Consider first the gas parcel in motion along a plate with no temperature 
gradient (Fig. 5a). The gas parcel starts at position 1 with a temperature T. The acoustic 


wave moves the parcel to the right a distance dx to position 2 and compresses it 


adiabatically, raising its temperature to T**+. The temperature of the gas is now higher 


than the local ternperature of the plate so heat Q flows from the gas to the plate, causing 


the plate to warn up at position 2 and the gas to cool to temperature T*. As the gas cools 
its volume decreases, so that work is done on the parcel by the surrounding gas. The gas 


parcel now moves back to position 1, undergoing adiabatic expansion. It arrives at 


position 1 with a temperature T-. which is lower than the local plate temperature. The 
result is a heat flow frorn the plate to the gas parcel. The increase in temperature of the 
gas parcel causes it to expand, resulting in work being done by the parcel on the 
surrounding gas. As this cycle is repeated, the constant dumping of heat at position 2 


and the removal of heat from position 1 results in a temperature gradient across the 


plate and heat transfer to the right along the plate. 








{a) Heat pump (b) Prime mover 


Figure 5; Fluid parcels in (a) a thermoacoustic heat pump and (b) a prime mover. 
Steps 1 and 3 are adiabats while 2 and 4 are isochors. 





Next, consider the cycle for a prime mover (Fig. 5b). Assume a temperature 
gradient (VT} is maintained along the plate by means of heat exchangers located at 
either end of the plate, with the hoi heat exchanger located to the right. Now when the 


gas parcel is compressed adiabatically and moved from position 1 to position 2, heat 


will flow from the plate to the gas parcel as long as T + dx(VT) > T*+t. This process will 


result in an expansion of the gas parcel as it is heated to a temperature T***, so that 


work is done on the surrounding gas. When the gas ts then adiabatically expanded and 


moved back to position 1 it will be at temperature T+. Heat must then flow from the gas 
to the plate and work will be done on the parcel by the surrounding gas. Heat is thus 
pumped from the hot heat exchanger to the cold heat exchanger. 

The actual process is sinusoidal rather than articulated as described above. 


As we Shall see later, the heat transfer takes place only within a thermal penetration 


depth of the surface of the plate. This thermal penetration depth (5x) is the distance heat 
can diffuse through the gas in a time equal to the reciprocal of the angular frequency at 
which the motion ts taking place. Thermoacoustic engines may be thought of as a long 
train of adjacent gas parcels, all within about a thermal penetration depth of the 
surface of the plate, that draw heat from the plate at one extreme of their motion, and 
deposit heat at the other extreme (Wheatley, Swift, and Migliori, 1986). 
2 The Single Plate 

The basic principles of thermoacoustic engines are inost Clearly illustrated 
by examining the simple case of a single, small solid plate aligned parallel to the 
direction of vibration of a standing acoustic wave. The presence of the plate in the 
standing wave modifies the wave, resulting in (1) a time-averaged heat flux near the 
surface of the plate, along the direction of the acoustic vibration and (2) the generation 
or absorption of real acouslic power near the surface of the plate. These two simple 


effects are the basis of all thermoacoustic engine phenomena (Swill, 1988). The 


il 





following discussion follows Swift's paper closely. The only significant difference ts 
that we have chosen to mcdel the acoustic pressure as a cosine function so that the 
resulting expression for fluid particle velocity contains a sine term. 

In the absence of the plate, the pressure and x-direction velocity of the 


standing wave have well known acoustic properties. Let the Muid pressure be given by 


p= pm + piel (5) 


where the acoustic pressure is given by 


Pl = Po cos(kx) =p} 5(x). {6) 


Figure 6: Geometry of the single plate thermoaccustic example. | 





The fluid particle velocity is fuund by applying the linear inviscid force equation 
(Kinsler and others, 1982) to Eq. 5. The result is 


u= ud, (7) 


where 


ul = -i(Po/pmc) sin(kx) =-iu15(x). (8) 


In the absence of the plate, the sound wave is adiabatic and has an 


oscillating temperature given by 


Tl = (Tmf/pmcp)p}. (9) 


For an ideal gas we know 


TmB/pmcp = (y-1)Tm/ypm. (10) 


Combining Eqs. 9 and 10 it can be shown that 


T1/Tm = I-1)/yip1/pm (11) 
for ideal gases. Thus the fractional temperature oscillations are of the same order of 
magnitude as the fractional pressure oscillations. 

The introduction of the plate into the standing wave modifies the 
unperturbed temperature oscillations due to the heat flow between the gas and the plate. 
The presence of the plale inodifies the temperature oscillations in both magnitude and 
phase, for fluid within about a thermal penetraticn depth of the plate surface. The 


thermal penetration depth in the fluid is given by 


8x = V(2xK/w). (12) 
At this point it is prudent to make several assumptions that will facilitate 
the calculations which follow. These are: 
(1) that steady state exists, 


(2) the plate is short compared to the acoustic wavelength and that it ts located far 
from both the pressure and velocity nodes so that both p1 and uj can be considered 


uniform over the plate, 


13 





(3) the fluid ts inviscid so that u) does not depend on y, 


(4) the plate has a large enough heat capacity per unit arca So that it's temperature 
does not change appreciably at the acoustic frequency, 

(5) the thermophysical properties of the plate material and the gas are 
ir tependent of temperature, 


(L, the plate has a mean temperature gradient in the x direction VTm and the plates 


thermal conductivity in the x direction is negligible, 

(7) the fluids thermal conductivity in the x direction is also negligible, 

(8) the mean fluid temperature is independent of y and is the same as that of the 
plate. 


Calculation of the oscillating fluid temperature is accomplished using only 
the first order terms in the general equation of heat transfer. expressing the entropy in 


terms of the pressure and temperature of the fluid as 


S1 = (cp/TmIT1 - (B/pm)p1. {13) 


and finally solving the resulting differential equation subject to the appropriate 


boundary conditions (Swift, 1988, p. 1151). The solution is 


Ti=[(TmB/pmcp)pis - (VTm/w)urS|{1 - eM +y/6x). (14) 


To understand the terms of this equation, first take the limit as 5x << y where the fluid 


makes negligible thermal contact with the plate. Under this condition 


Tl {TmB/pmcp)p1§ - (VTm/)u15. (15) 


The first term is due to the adiabatic compression and expansion of the gas. The second 


term is a result of the mean temperature gradient in the fluid and the fact that the fluid 


oscillates in the x direction with displacement amplitude ui%/w. Setting T1 = 0 and 


solving Eq. 15 for VTm. we see that there is a critical mean temperature gradient 


VTert = TmBwp1$/pincpais (16) 


for which the temperature oscillations at a point are zero. This critical temperature 


gradient is important because !t marks the boundary between the heat pump and the 





prime mover regions of operation. Further manipulation of Eq. 16, using Eqs. 6 and 8, 


and the generalization of Eq. 10 to an arbitrary fluid leads to the useful result 


VTcrit = [(y-1)/TmB](Tmk) cot(kx) (17) 


The exponential term in Eq. 14 is complex. It approaches zero for y << 5x 


where the plate imposes the boundary condition T} = 0. When y = 5x it has a substantial 
imaginary part. This phase shift in the oscillating temperature of the standing wave 


due to the presence of the plate leads directly to the time averaged heat flow in the x 


direction, 


Q Heat Flux 
Since we have assumed that thermal conductivity in both the plate and 
the gas is negligible in the x direction, the only way that heat can be transported in the x 


direction is by the hydrodynamic transport of entropy, thus the heat flux per area is 


given by 

q2 = Tmpm<s]ul >. (18) 
Substituting Eq. 13 for s] and simplifying results in ar expression for heat flow as a 
function of y. The total heat flux Q2 along the plate in the x direction can then be found 


by integrating q2 over the y-z plane. The resulting expression is 


Q2 = -(1/4)16xTmBp1Su1S(r-1) (19) 


where T = VTm/VTenit is the ratio of actual temperature gradient to critical temperature 


gradient. The term Tmf ts nearly 1 for ideal gases and the term [15x is the area which is 


thermodynamically active in a plane perpendicular to the plate. The heat flux will be 


zero if the plate is positioned at either a pressure or a velocity node of the standing wave 


which is readily understood from the product p15u15. Finally we see that if VIm= VTerit 
then [ = 1 and no heat flux occures. If VTm < VTcrit . then r—1 < 0 and the heat flux is 


toward the velocity node. If VI'm > VTcrit , then T-1 > O and the heat fNux ts toward the 
pressure node. If suitable heat exchangers are positioned at the ends of the plate, the 
heat flux carries heat from one heat exchanger to the other. 
b. Work Flux 
The work flux. or acoustic power is derived in a similar manner from 
thermodynamic principles. The work done by a differential volume element of Muid in 


going from a volume dx dy dz to a volume dx dy dz + dV is 


dw = p dV = - (p/p) dx dy dz dp, (20) 


so the power per unit volume is 


w = - (p/p) dp/dt. (21) 
Since, for oscillatory motion, we are interested in the time average work flux we use Eq. 


5 for p and 


dp/dt = twp) + u1{dpm/ax) 


w2 = - (wpm) <ipip1> (23) 


for the average acoustic power produced (or absorbed) per unit volume. After expressing 


pl in terms of T) and p) and substituting intc Eq. 22 we finally arrive at the result 


w2 = - (1/2) oBpi§ Im{T)). (24) 
Just as with the heat flux per area, only the imaginary part of the temperature 
oscillations contributes to the acoustic power produced per unit volume. The fluid 


within about a thermal penetration depth of the surface of the plate “breaths” because of 





thermal expansion or contraction, with the correct time phase with respect to the 
oscillating pressure to do or absorb net work. 
By integrating w2 over all space, using Eq. 14 for T), the total acoustic power 


produced is found to be 


W2 = (1/4)118«Ax(TmB*w/pmcep)(p18)2(r-2). (25) 


This acoustic power is proportional to the volume of the fluid, []5xAx, that is within a 
thermal penetration depth of the plate. It is also quadratic in the acoustic amplitude 


and vanishes at the pressure nodes. Finally it is proportional to the same temperature 


gradient factor T-1 as appeared in the expression for the total heat flux. Whether this 
power increases the amplitude of the standing wave, is radiated away to infinity, or 
flows through an acoustic-to-electric transducer to generate electric power depends on 


the details of the resonator. 
c« Efficiency 


The efficiency of the single plate engine, used as a prime mover {s given 
by 1 = W/QH. For a short plale, Q2 is essentially constant with x, so QH = Qc = Q2. 


Substituting for W2 and Q2 from Eqs. 19 and 25 above results in the expression 


n= AxBwp1$/pmepulS. (26) 


Using the definition of VTcrit, this equation can be recast Into the form 


He 1Ke/T. (27) 


Thus we see thal as T - 1, the single plate engine approaches Carnot's efficiency, but 
this is also the point where the vower output of the engine approaches zero. 
3 =. The Short Engine 
.. the previous sect‘on, thermoacoustic effects were examined for the simple 


case of 7s rdle plate in an acoustic standing wave. In order to examine these effects 
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without a great deal of distractions, several simplifying assumptions were made. The 
central thermoacoustic effects that were revealed are: 
(1) The thermal boundary condition imposed on the fluid by the plate causes a 
phase shift of the oscillating temperature of the fluid within a thermal 
penetration depth of the plate. 


<) Heat flows hydrodynamically, parallel to the plate in the fluid about a thermal 
penetration depth from the plate. 


(3) Acoustic power is absorbed or produced by the fluid about a thermal 
penetration depth from the plate. 

(4) The direction of heat flow. and whether power is produced or absorbed, depends 
on the magnitude of the mean-iemperature gradient relative to the critical 
temperature gradient. 

We now turn our attention to the more realistic model of an engine built with 
a stack of parallel plates. The two sets of coordinate axis used in the detailed 
calculations (Swift, 1988. pp.1176-1179) which are summarized below are indicated in 
Fig. 7. We shall also discard many of our previous simplifying assumptions although we 
Shall retain the assumptions that the temperature spanned is much less than the 
absolute temperature and that the length of the plates is much less than an acoustic 
wavclength, 

Swift's derivation proceeds as follows. Consider a parallel plate engine 
geometry with the acoustic vibration along the x axis and the y axis perpendicular to 
the planes of the parallel plates. The equation of motion for the Nuid is 

plav/at + (v-V)v) 
= -Vp + pV2v + (&+/3)V(V>v), (28) 
with the boundary condition v = O at the surface of the plates. The continuity equation 


for the fluid is 


dp/at + V:(pv) = 0. (29) 


The equations for heat flow in the fluid and solid are 


pTids/dt+v° Vs) = V°(KVT) 


+ (terms quadratic in velocity) (30) 
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and PsCsdTs/dt = KsV2Ts, (31) 


with boundary conditions T=Ts and KdT/dy=KsdTs/dy at the fluid-solid interface. The 
above equations are linearized, assuming that each of the variables oscillates at 


angular frequency w: 


plates 


Figure 7: Geometry used for multiplate engine calculations. (a) Overview of short 
engine in resonator, and (b) expanded view of plates. Each plate has thickness 2L 
and is separated from adjacent plates by a distance 2yo. 
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P= DPm+ pixel, (32) 


P= Pmix) + pilx.y)el (33) 
¥ =lurx.yles ur (xye!, (34) 
T = Tmlx) + Ti(x.y)!o, (35) 
Ts = Tm(x) + Ts} (x.y), (36) 
S = Sm(x) + Si(x,y)!@t, (37) 


Integrating the resulting equations with respect to y yields a wave equation 


for pi(x) in terms of dTm/dx and material properties and geometry, and an equation for 


energy flux H2 (which includes both heat and work) along x in terms of p1(x), Tm(x), 


material properties, and geometry. The wave equation is 


{1+{y-1)fk/(1+es))p1 + (pme2/2)d/dx{{(1-fv)/pmidp1 /dx} 


- (Bo2/w?){(fe-fv) /{(1-9)(1 +€5)}(dTm/dx)(dp1/dx) = 0, (38) 
where 
fy = tanh{(1+i)yo/8vI/{(1+i)yo/ Sv). (39) 
fx = tanh{{1+i)yo/5x)/{(1 +i)yo/ 8x}. (40) 
Es = pmcpdxtanh|{(1+iyo/8xl/pscsSstanh[(1+1)L/6s], (41) 
by = V(2v/a). (42) 
Sx = V(2K/@), (43) 
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and ds = V(2xs/u). (44) 
The term 8s is the thermal penetration depth in the solid, and is analogous to the 
thermal penetration depth in the fluid discussed previously, while 5y is the viscous 


penetration depth of the fluid. The viscous penetration depth is roughly the distance 


shear diffuses in a time 1/w. The equation for energy flux is 
H2 = (Myo/2upm)Im{(dp1°/dx)pil1-fv"’-Tm Blfe-fv")}/{(1+0)(1+es)} I] 
+ (Myocp/2u3pm(1-9))(dTm/dx)(dp1 /dx)(dp1" /dx) 


X Imffe' + (ff) +eslfy /fx) /{(1+0)(1+es)}] 


- TiyoK+LKs}(dTm/ dx). (45) 

In order to make contact with the single plate resulis, we make three 
assumptions designed to simplify the equations. These assumptions are zero viscosity. 
the “boundary-layer” approximation, and the “short-stack" approximation. (The 
assumpton that the fluid viscosity is zero is not a good approximation in most realistic 
situations, requiring that the viscous terms bi ut back into the equations at a later 


time.) With the fluid viscosity set equal to zero the terms o and fy both vanish. 


The “boundary-layer” approximation consists of the assumptions that 


yo>>&x and that L>>8s. The goal of these assumptions {ts to set the hyperbolic tangents 


equal to unity. This is reasonable for cases where the value of the dimension of interest 


is at least twice the value of the associated boundary layers (i.e. yo 2 28« )[Swift, 1988). 


The “short-stack” approximation assumes that the stack is short enough 


that it does not appreciably disturb the standing wave, so that we can take 


Pl = Po coslkx) = p15 (46) 
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ul = -i(1+L/yo)(Po/pmc)sin(ky) = -iu1S. (47) 
These are the same expressions used in the single plate calculations except for the term 


1+L/yo which is required for continuity of volume velocity. It is further assumed that 


the stack is short enough that p1$ and u15 may be regarded as independent of x within 
the stack, and that AT=TH-TC <<Tm, so that thermophysical properties can be assumed 


independent of Tn and hence of x. 


With these assumptons, and using the equation 


dpi/dx = -wpmu}$, (48) 


(Swift, 1988, p. 1158) we obtain 
H2 = -(1/4)115xlTmBpiSui§/(1+es)][F-1) 
- MlyoK+LKs)(dTm/ dx) (49) 
for energy flux, and 
P1 + (pme2/w*)d/dx{(1/pm)dpi /dx} 


= [(y-18xp1/((1+)(1+eslyol|(T - 1) (50) 
for the wave equation. 


The acoustic power generated (cr absorbed) in the stack may be determined 


by multiplying the change in average acoustic intensity between the two ends of the 


stack by the fluid cross section at the stack exit. The resulting equation is 
W2~= = Tiyol<p1ul>right - <p)ul>left] 
= [yoAxd(<pju}>)/dx. 


Marupulaling Eq. 51. using Eqs. 48 and 50, we find that 





W2 = (1/4)T18xAxl(y- 1)@/{pme(1 +es)}](p18)2(T- 1) (52) 
We pause now to make contact with the simple results of the single plate 


calculations. We begin by examining the energy flux equation. The first term on the 


right in Eq. 49 differs from the right side of Eq. 19 by a factor 1/(1+es) which accounts 


for the finite heat capacity of the plates and modifies the boundary condition for 
temperature at the fuid-solid interface. The second term on the right side of Eq. 49 
accounts for ordinary heat conduction in the plates and the fluid. Upon examination of 


the equations for work flux (Eqs. 25 and 52) we find that the difference is again the 
factor 1/(1+€s) described above. 


We now include viscosity and repeat the above calculations, but now use the 


relationship 

api/dx = -wpmu1S/(1-fy) (53) 
(Swift, 1988, p. 1159). Notice that dpi/dx is now complex due to the presence of viscosity 
in the fluid. (Note that Eqs. 46 and 47 may still be used for this case as long as we now let 


u}% represent the mean fluid velocity in the x direction. Due to the presence of viscosity. 


the x component of the fluid velocity is a function of y, going from zero at the surface of 
the plates to some maximum value at a distance on the order of a viscous penetration 


depth away from the plates.) We now obtain 


H2 = -1/41] &klTmBp1$u1$/{(1+0)(1 +€s)(1-5v/yotbv2 /2yo2)}] 


X [P1+Vo+0+0es)/(1+Va) - (1+V0-5y/yo)] 


- T(yoK+LKs)(d fm /dx) (54) 


for energy flux, arid 





Pl + (epme2/w2)d/dx{((1-f)/pmidpi /dx) 


= (7-DSkp1/((1+i)(1 +es)yo}[P/i(1+Vo)(1-fv)} - 1 (55) 
for the wave equation. 


The work flux is calculated in the same manner as before and the result ts 
W2 = (1/4)TT5xAxi(y- 1)o(p18)2 /(ome(1+es)}] 
X (P/{(1+V0)(1-5v/yor8v2 /2yo%)} - 1] 


- (1/2) [TSvaxlopm(u1s)2/(1-8v/yorbv2 /2yo2). (56) 
It is clear that the inclusion of viscosity in the calculations adds significant complexity 


to the results, and the intuitive understanding of the components of the equations is 


lost. 





Il. DESCRIPTION O * APPARATUS 


The apparatus which we are modeling in this thesis is a thermoacoustic prime 
Mover, enclosed in a rigid tube resonator, driving an acoustic-to-electric power 
transducer. Figure 8 is a general schematic of the device. Specific values of the 
dimensions of the prime mover are contained in Table 1. 

The transducer proposed for use in this application is a linear alternator. Design 


requirements submitted for bid by Dr. Greg Swift of L.A.N.L. include 1 kW of power 


a) 


Thermoacoustic 
stack assembl Resonator 


Transducer 


Hot heat Cold heat 
exchanger Stack = exchanger 


Figure 8: Schmatic representation of a thermoacoustic power generator. (a) 
Overview of system. (b) Details of stack assembly. 
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| Hot heat Cold heat 
exchanger exchanger 


Characteristic 


Inlet position 
relative to rigid 
end (m) 


Exit position 
strike to rigid 





Table 1: Summary of geometry of thermoacoustic stack assembly an 

important thermodynamic parameters. 
delivered to an electric load resistance of 3 to 100 ohms with at least 60% efficiency 
when driven by a 4 bar peak-to-peak pressure amplitude at 120 Hz in helium at a mean 
pressure of 20 bar. The details of the proposals received in response to the request for 
quotation contain proprietary information. The transducer model used in this thesis 
was designed with the assistance of Nick Vitale of MTI. Although it does not represent 
the actual design proposal submitted by MTI, it should provide reasonuble 
approximations for the performance of the prime mover-transducer system. The 


transducer may be modeled by the equivalent circuit indicated in Fig. 9. The moving 


mass of the transducer is 0.766 Kg, the spring constant is 1.43E5 Kg/s2, the magnetic 
field strength is 0.97 Tesla, and the coil has a total effective length of 32.6 m. The 
mechanical resistance of the transducer ts 2.0 Kg/s. Although the actual transducer will 
have some physicul limitation to its displacement, we have not imposed any such 


limitation on our model. 








Transducer 






where: Ri =562 ohms 
Fe =0844 ohms 
L=240mH 
C=732pF 

Road = variable 













Figure 9: Equivalent circuit model of linear alternator acoustic-to-electric power 
transducer. 





IV. NUMERICAL MODEL 


In order to formulate a program capable of modeling the performance of our 
thermoacoustic generator, we must identify the parameters which may be allowed to 
vary during the operation of the device, and then develop a system of equations that 
will relate these variables. The geometry of the device will be fixed. It is assumed that 
the thermodynamic properties of the material used in the constriction of the device 
will not vary significantly over the range of pressures and temperatures at which the 
device will operate. It is further assumed that the temperature of the coid heat 
exchanger will be constant. The working fluid used will be a mixture of helium and 
argon gases. The mixture of these gases will be determined so as to achieve a 120 Hz 
operating frequency at a median load resistance of 13 ohms. The parameters that will 
be allowed to vary in the model are the acoustic pressure amplitude. angular frequency, 
and the temperature of the hot heat exchanger, for a given input heat flux and load 
resistance. Variation of the hot end temperature will result in variation of the 
thermodynamic properties of the working fluid and of the sound speed in the fluid. 

Before continuing, one additional assumption will be made concerning the 
performance of the system; specifically, that no attenuation occurs in the resonator. 
All acoustic power generated in the stack must therefore be removed at the transducer. 
This assumption is made to simplify the calcuiations, but is not truely valid since a 
significant fraction of the acoustic power would be attenuated due to thermal and 
viscous effects at the walls of the resonator. 

We now begin to develop the equations which will determine the performance of 
the system using the short stack and boundary layer approximations discussed 
previously. Starting at the rigid end of the resonator, and assuming a standing wave, we 


know that the pressure amplitude must be at a maximum and the fluid particle velocity 
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must be zero. We can therefore describe the pressure and fluid particle velocity in the 
region between the rigid end and the stack assembly using Eqs. 6 and 8, which are 


repeated below for ease of reference. 


p(x) = Po cos{kx) (6) 


u(x) = -i(Po/pmc) sin(kx) (8) 
The short stack approximation allows us to assume that the pressure amplitude in the 


stack is constant. Using the mean Stack position to determine the value of this term we 


find 


plxm) = p15 = Po cos(Kkmxm) (57) 
where km = w/cm. (58) 


and cm is the mean sound speed in the system. The mean fluid particle velocity at the 


mean stack posilion is given by 


u(Xm) = -iu1 = -ig(Po/pmc) sin(kmxm), (59) 
where g = 1+L/yo is the geometry factor necessary to maintain continuity of volume 


velocity. It must be noted that g has different values for the stack and the two heat 
exchangers. 
The change in acoustic pressure due to the stack and heat exchangers is obtained 


from Eq. 53. The resulting equations are 


dpht = -Powdxhtght sin(kmxm)/{em(1-fv,ht)}. (60} 
dpst = -Powdxstgst siin{kmxm)/(cm(1-fv,st)}. (61) 
and dped = -Powdxedged sin(kmxm)/{cm(1-fv.ed)}. (62) 
where fv = 1/{(1 +i)yo/&). (63) 


Nolte the subscripts indicating the parameters which have different values depending 
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on the location in the stack assembly. The acoustic pressure at the exit of the stack 
assembly is then the sum of the pressure changes in the stack and the pressure at the 


stack entrance: 
plx4) = p(x1)+dpht+dpst+dped (64) 
where p(x1) = Po cos(khtx1). (65) 
The change in mean fluid particle velocity is obtained from the wave equation 


(Eq. 55) by substituting for dp1/dx from Eq. 53 and solving for du. The resulting 


equation is 
du1=ilop1dx/(pmc2)] 


X [Cy Sx /((1 +) +es)yo) (IP / (C1 +¥o)(1-fv))-1)- 1). (66) 


The changes in mean fluid particle velocity in the heat exchangers and the stack are 


then 
duht=ilwPodxht cos(Kmxm) /(pmem2)] 


X [{(y- Sx/((1 +i) 1+es, htyo. hor /((1+vo)(1-fy,nv)-1)-2], (67) 


dust=ilPodxst cos(kmxm) /(pmcm2]] 


X [l(y- 18x /((1 +1) +s, s0yo.st}(P/((+¥0)(1-fv,st))-1]-1). (68) 


and 
ducd=![wPodxed cos(kmxm)/(pmem2}} 


X Illy 18x /((1+1)(1+es,cd)yo.cd)}[P/((1 +vo)(1 -fy,cd))-1)-1). (69) 


The fluid particle velocity at the stack exit is 


u(x4) = u(x})+dunt/ght+dust/gst+ducd/god (70) 


where u(x)) = -i(Po/tmcht) sin(khtx1). (71) 








The geometry factors are required for continuity of volume velocity as before. 
Assume now that the acoustic wave in the region between the stack assembly and 


the transducer has the form of left and right traveling plane waves described by 


p{x,t) = Aellot-lod 4 Betlat+kag (72) 
where A and B are complex coefficients given by 

A=a+ia (73) 
and B=b+9p. (74) 


The associated fluid particle velocity is found from Euler's equation to be 


u(x,t) = (A/pmecdell@-d ~ (B/pmecalellot+h0 (75) 


Dropping the time dependence of Eqs. 72 and 75, and expanding in terms of sines and 
cosines yields equations for acoustic pressure and fluid particle velocity in the 


resonator between the stack assembly and the transducer of the form 


p(x) = (a+b)cos(Kedx)+(a—-B)sin(Kedx) 


+ i[(a+B)cos(kedx)-(a-b)sin(kedx)] (76) 


u(x) = (1/pmcecd){(a-b)cos(kedx)+(a+B)sin(kedx) 


+ i[(a--P)cos(kedx)-(a+b)sin(kKedx)}) (77) 


One set of boundary conditions for Eqs. 76 and 77 are given by our expressions for 
pressure and velocity at the stack assembly exit. This requires expansion of the 


expressions into real and tmaginary parts so that we may form the equalities 


(a+b)cos(Kedx)+(a-B)sin(kedx) = Re{p(xi)+dpht+dpst+dpcd}, (78) 
(a+B)cos(kedx)-(a-b)sin(kcedx) = Im(p(x1)+dpht+dpst+dpcd). (79) 


(1 /pmced)}{(a-b)cos(Kedx)+ (a+$)sin(kedx)} 
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= Re(u(x))+duht/ght+dust/gst+ducd/£cd}. (80) 
and (1/pmCcd){(a~B)cos(Kedx)-(a+b)sin(kedx)} 


= Im(ulx})+duht/gnt+dust/gst+ducd/ged}. (81) 


The boundary condition at the transducer required for resonance is that the 


reactance of the combined system must vanish: 
Im{Zr + Zm) = O. (82) 


The mechanical impedance of the acoustic wave at the transducer position is given by 


Zr = (nd2/4)[plxy/ulxs)]. (83) 
The mechanical impedance of the transducer is found by applying Newton's second law 
to the transducer. The force acting on the transducer by the acoustic pressure ts 

Ft = (nd2/4)p(xi) (84a) 
Applying Newton's second law to the transducer gives 

Fi =mia+s$+Rmu +BLel. (84b) 


This equation may be re: ast by recognizing that 


a= 1wu, (85) 
0 = -iu/o, (86) 
and 1= Bleu/Ze, (87) 
where Ze = 11/Ri+1/{R2+Rload+ilwL-1/aC)})"}. (88) 


The result is that the mechanical impedance of the transducer may be written as 


Zm = Ft/u = Rm +(BLc)2/Ze +i(wmt-s/w). (89) 


We may find two more equations, necessary to fully constrain the 7 variables of 


interest. by recognizing that the conditions in the system are directly effected by the 


amount of heat flux in the stack, and that the power propagating from the stack 





assembly must be removed at the transducer during steady state operation. Our sixth 


equation is therefore Eq. 54. The seventh equation is 


Tac - Im - Me = 0 (90) 

where Tac = (xnd2,/8)Re(pu’), (91) 7 
Tm = (1/2)Rm{uu‘}, (92) 

and Tle = (1/2)Re{I*V). (93) 


The voltage produced by the transducer is 
V =Bleu. (94) 


Before describing the method used to firid solutions to these equations, it is useful 


to stummarize them and the variables of interest. The seven equations are 
(1)  (a+bjcos(kedx)+(a-B)sin(Kedx) = Re(p(x1)+dpht+dpst+dpcd}, 

(2) (a+B)cus(kedx)-(a-b)sin(kedx) = Im{p(x1)+dphi+dpst+dpcd). 

(3) (1/pmced){(a-b)cos(Kedx)+(u+B)sin(Kedx)} 


= Reju(x})+duht/ght+dust/gst+ducd/ficd), 


(4) (1/pmccd){(a-B)cos(Kedx)-(a+b)sin(kKedx)} 


= Im{u(x1)+duht/ght+dust/gs1+ducd/gcd). 


(5) Im{Zr+Zm}=0 , 


(6) He = -1/4718«(TmBp1$ui§/{(1+0)(1 +s 1-5v/yorby2/2yo))] 


X [P[{1+vo+a+0€s}/((1+va)} - (14V0-8v/yo)] 


- MlyoK+LKs)(dTm/dx), 





and 

(7) lac - Tm- We = 0. 

The variables for which we are to solve these equations simultaneously are the five 
pressure coefficients Po, a, a, b, and B, the angular frequency , and the temperature of 


the hot heat exchanger Tht. In order to accomplish this task effectively, it is necessary 


to expand the above equations in terms of the parameters of interest and to solve 
analytically for the real and imaginary parts of the equations indicated, This was 
accomplished using an algebraic manipulation software system called Derive: a 
Mathematical Assistant for Personal Computers. The details and results of the 
calculations are to long to present here. 

Before we begin to solve these equations we must know the values of all the 
parameters that are not variables. The parameters that describe the geometry of the 
system and the characteristics of the load will be fixed at the beginning of the process. 
The parameters that describe the working fluid are dependent on the gas mix used and 
on the temperature at the hot heat exchanger. This temperature effects the values of the 
viscosity, thermal conductivity, and Prandt] number of the gas mix as well as the 
sound speed in the system. which in turn effects the gas mix that will be chosen in order 
to have resonance at the desired frequency. It ts necessary therefore, to begin the 
calculations by guessing a value of the hot end temperature and, based on this guess, 
calculate the required gas parameters. We must then solve the equations to find a new 
hot end temperature, and then adjust the values of the gas parameters in an iterative 


cycle, until a self-consistent solution ts obtained. 


We begin the process by “guessing” a value for Tht. Using this value we calculate the 


mean temperature gradient VTm, the mean temperature in the stack Tm, the average 


system temperature Tavg, and the crilical temperature gradient VTcrit using the 





relations 


VTm = (Tht-Ted) /dxst. (95) 
Tm = (Tht+Tcd)/2, (96) 
Tavg = (x2Tht +dxstTm +(x5-x3}Ted)/x5. (97) 
and VTcrit = {(y-1)/TmB](mkm) cot(kmxm), (17) 
where k = w/cm = 2/Xxt. (98) 


For this first approximation it is assurned that TmB = ) and y=1.67. Using the result of 


the above calculations we then determine the normalized temperature gradient 


r=VTm/VTent. 

With the temperature parameters calculated above we now determine the first 
estimate of the gas mix required to obtain resonance at the desired frequency fo. This is 
accompUshed by using the expression fc. sound speed in a perfect gas 

c = V(yrT) (99) 
and estimating the required sound speed as c = foA = 2foxt. The gas constant is given by 
r= R/(FheWhe + FarWar) (100) 


where R is the universal gas constant, F is the mole fraction of the gas indicated by 
subscript that is present, and W is the molar weight of the gas indicated by subscript. 


Since only the two gases are present in the system, we may write 
Far = 1 - Fhe. (101) 
Substituting Eqs. 100 and 10) into Eq. 99. and using c = 2foxt, we solve for the mole 


fraction of helium required: 


Fhe = 1/(Whe-War)[yRTavg/(2foxt)? - War). (102) 





From these results we calculate the molar weight of the binary gas mixture for future 
use: 
Wmix = FheWhe + FarWar. (103) 

Having established an approximation for the required gas mix, we next determine 
the isobaric and isochoric specific heat capacities of the gases at Tavg. This is done by 
two subroutines called ‘HELIUM’ and 'ARGON'. These subroutines are modifications of 
programs, written by Robert D. Mccarty at the National Bureau of Standards, to 
calculate the thermodynamic properties of these and other gases, using a form of the 
virial equation of state. The modifications to the programs consisted primarily of 
eliminating functions not required for this application, such as output to files and 
calculation of additional thermodynamic properties (internal energy, entropy, 
enthalpy. etc..), and modifying the input requirements so that the subroutines read 
temperature and pressure (the input parameters required) from the main program 
instead of from the screen in an interactive mode. In addition, the upper temperature 
limit in ‘ARGON’ was raised from 400 K to 1000 K. Comparison of data calculated in this 
extended temperature region was found to be in agreement with tabulated data 
contained in available tables of thermodynamic properties for argon. Using the output 
from the subroutines, the isobaric and isochoric specific heat capacities of the binary 


gas mix are determined: 

Cp.mix = (Cp,heFheWhe + Cp,arFarWar)/Wmix (104) 
and Cy.mux = (Cv.heFheWhe + Cv,arFarWar)/Wmix. (105) 
From these, we then calculate a new value of the ratio of heat capacities, y = 


Cp.mix/Cv,mix. If y varies significantly from the initial value of 1.67, the gas mix is 


calculated again, using this new value, until a self consistent value is achieved. 


Now that a reasonable gas mix has been established, and a value of the ratio of 


heat capacities determined, we can calculate the mean gas density and acoustic sound 





speed in the resonator, assuming ideal gas behavior. The mean gas density is found 


using the ideal gas relation 


Pm = pmWmix/RTavg. (106) 
The acoustic sound speed in the resonator is found using Eq. 99. 

Of more interest are the gas properties in the stack assembly, where the 
thermoacoustic effects are occurring. To determine these, we again use the subroutines 
‘HELIUM’ and 'ARGON’, but with Tm as the input temperature. We obtain the specific 
heat capacities of the two gases, calculate the specific heat capacities of the binary 
mixture, and determine y and the acoustic sound speed for the stack assembly as 
described above for the resonator. To determine the dynamic viscosity, thermal 
conductivity, and Prandtl number of the binary gas mixture, we call subroutine 
‘PRANTL’. PRANTL is a modified version of a program written by a student at the Naval 
Postgraduate School (Susalla, 1988) to calculate these same parameters for 
thermoacoustic refrigerator design calculations. The original program calculated the 
dynamic viscosity, thermal conductivity, and Prandtl number for helium-argon and 
helium-xenon gas mixes at only two temperatures. The program was modified to allow 
for calculations over a range of temperatures, but only for helium-argon mixtures. The 
input parameters required for the subroutine are the temperature of interest and the 
{sobaric specific heat capacities of the two gases. 

Having now obtained the necessary parameters based on our ‘irst estimate of the 
hot end temperature, we are nearly ready to solve our seven equations. The seven 
equations which we must solve simultaneously are non-linear in the variables of 
interest, and can not easily be solved analytically. We chose instead, to use a readily 
available Fortran Library subroutine to find the solution for us. The subroutine used is 
the IMSL Library subroutine NEQNF. NEQNF is a single precision subroutine which 
uses The Levenberg-Marquardt algorithm and a finite-difference approximation to the 


Jacobian to solve a system of nonlinear equations. 
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In order to find the solution, the subroutine requires initial "guesses" for the 
values of the variables. We have established a "guess" for the hot heat exchanger 


temperature in order to estimate the gas properties above. The "guess" for the angular 


frequency is obviously = 2nfo. What is lacking are initial values of the pressure 
coefficients. 


In previous work with the 5 inch engine at L.AN.L., operating it in a simple rigid- 


rigid resonator, acoustic pressures on the order of 10% of the mean pressure have been 


observed. We will take this value for Po. Values for a, a, b, and B may be obtained by 


analytically solving equations similar to Eqs. 78-81 derived from the inviscid 
relationships for the short engine. The resulting equations are too long to be repeated 
here. Initial "guess” values for a, a, b, and § are calculated in subroutine ‘FIRST. 

Having determined all of the required input parameters, we now call the library 
subroutine to solve our system of equations for the seven “unknown” variables. If the 
resulting new hot end temperature differs significantly from the initial "guess" it is 
necessary to redo the calculations of the temperature dependent parameters and gas 
parameters (excluding the gas mix) and solve again. This process is repeated until the 
values of the hot heat exchanger temperature do not change significantly from one 


iteration to the next. 
Once a consistent value of Tht is determined, we compare the resulting angular 


frequency to our desired operating frequency. If the calculated frequency is higher than 
the design point, we increase the argon fraction, thereby reducing the acoustic sound 
speed, and repeat the {iteration process. Conversly, if the frequency is too low, the argon 
fraction is reduced, and the iteration process repeated. These corrections are repeated 
until a gas mix is found which puts use near the desired operating frequency. 


Finally we compare the electric power generated by the transducer (not that 


delivered to the load) with a desired maximum value. If the power generated differs 





from the design value, H2 is adjusted and the entire iteration process is repeated. 


The ultimate result of the above process is a gas mix that will result in operation i 
of the systern at the desired frequency, and provide the desired power output, for a given 
load resistance and a particular heat flux input. 

It should be noted here that the above iterative processes are actually carried out 
in two stages. The first stage. accomplished in subroutine ‘INITIAL’, provides a rough 
estimate of the required values. Subroutine ‘REFINE’ takes these rough estimates and 
refines them to the final values. 

Once the gas mix is determined, it is fixed and the system of equations may be 
solved repeatedly for a variety of input heat fluxes and load resistances to determine 
the response of the system about the design point. The instantaneous power delivered to 


the load by the transducer is given by 

Th = V1. (107) 
The current flowing through the load resistor is given by 

I = (V/Z) sin(wt+9) (108) 
(Ohanian, 1985) where V is the voltage produced by the transducer given by Eq. 137, 9 is 
the phase difference between the voltage and the current given by 

$ = tan" }{(1/wC-wL)/(R2 + Rload)]. (109) 
and Z is the impedance of the branch containing the load resistor given by the relation 

Z = R2 + Rioad -i(@L-1/0C). (110) 
After appropriate substitution and reduction we find that the time average power 
delivered to the load is 

Toad = 1/2(BLc)*Rioad{u*u/Z*Z) cos(9). (111) 
The term cos() is very nearly unity for the values of L, C, and w of interest, and is 


neglected in the calculations of [load in the program. We may now calculate the total 
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system efficiency as 


1 = Ilload/H2. (112) 





V. RESULTS 


In this chapter we will present the results of calculations made for various load 
resistances, input heat fluxes, and transducer positions. The trends of some important 
thermoacoustic variables are presented graphically. Since no experimental device 
comparable to the one modeled has been built, no comparison of theoretical and 
experimental data is possible. The complexity of the governing relationships makes 


Intuitive prediction of the trends difficult. 


A. VARIATION IN LOAD RESISTANCE AND INPUT HEAT FLUX 

In order to investigate the effect of varying the load resistance and input heat 
flux, we first fix the transducer position at 4.0 m. We then calculate the gas mix 
required for resonance at 120 Hz. with a nominal load resistance of 13.0 Q. The result 
of these calculations is a gas mix of 97% helium and 3% argon. 


Table 2 contains some representative values of important parameters calculated 


Load Angular Hot Useful Overall Acoustic 

Resistance Frequency End Power Efficiency Pressure 

(Q) (rad/s) Temp. Generated {%) Amplitude 
(Bar) 





Table 2: Values of critical parameters as a function of load resistance for a 97% 
helium gas mix. The transducer is at 4.0 m and the input heat flux is constant at 5 
. kw. 
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for an input heat flux of 5 kW. Figures 10 through 14 indicate graphically the general 
trends of these parameters as the load is increased. Figure 15 shows the variation in 
acoustic pressure amplitude as a function of input heat flux for a number of load 
resistances. The output of the program indicated no variation in frequency, hot end 
temperature, or efficiency as the iuput heat flux was changed © ‘a given load resistance 


so no plots of these relationships are presented. 


B VARIATION IN TRANSDUCER POSITION 

Calculations were also done for transducer positions of 3.5, 3.0, and 2.5 meters. 
Table 3 and the figures which follow indicate the general trends in the parameters of 
interest as the transducer is moved closer to the Uermoacoustic stack. 

We clearly expect the argon percentage to increase so that we maintain a resonant 
system. It is unclear what effect this will have on the overall performance of the 
system. Changing the ratio of argon to helium will certainly change the 
thermodynamic properties of the gas mix, but whether that change will be adverse or 
beneficial is not intuitively obvious. 

Examining the following figures, we see that the normalized temperature 
gradient increases, the pressure amplitude increases, and the system efficiency 
decreases as the resonator is shortened for the same stack geometry. 

The increase in the normalized temperature gradient as the transducer is moved 


closer to the stack assembly may be understood by noting that the concurrent change in 


Helium Thermal Viscous Acoustic Norm. 

Fraction Penetration Penetration Pressure Temp. 

(m) (mole %) Depth Depth Amplitude Grad. 
(mm) 
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the gas mix results in a decrease in the sound speed and thermal penetration depth of 
the gas. These decreases result in an increase ii the temperature gradient in order to 


ro produce the same amount of acoustic power. 


en A mt He penn 
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Angular Frequency vs. Load Resistance 


Ang. 
Freq. 
(rad/s) 


10 $1 12 13 14 15 16 17 
Load Resistance (ohms) 


Figure 10: Plot of the variation of angular frequency with increasing load 
resistance for a 97% helium gas mix. The transducer is at 4.0 m and the input heat 
O kV 


Hot End Temp. vs. Load Resistance 
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Figure 11: Variation of hot heat exchanger temperature with load resistance for Xt 
= 4.0m., H2 = 5kwW. 





Power Delivered vs. Load Resistance 
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Figure 12: Piot of power dissipated in the load resistor as a function of the load 
resistance at a constant input heat flux of 5kW. 


System Efficiency vs. Load Resistance 
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Figure 13: Variation in system efficiency as a function of load resistance at 
constant heat flux.: 











Acoustic Press. Amp. vs. Load Resistance 
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14: Plot of acoustic pressure amplitude as a function of load resistance. | 








00 Pressure amplitude vs. Heat Flux 






ooree aa ae ew wt sw DD 






2500 3500 4500 5500 6500 7500 8500 9500 
input Heat Fiux 
@Re13 wR<16 







oR=8 o R=10 





Figure 15: Acoustic pressure amplitude as a function of input heat flux for various 
load resistances. 





Helium Fraction vs. Xdcr Position 
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Figure 16: Mole fraction of helium required for various transducer positions for 
13 ohm load resistance. 








230 Penetration Depths vs. Helium Fraction 
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Figure 17: Variation of penetration depths with helium fraction. 
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Acoustic Press. Amp. vs. Xdcr Position 
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Figure 18: Variation in acoustic pressure amplitude as a function of transducer 
osition to produce 1.1 kW of total electric power. 








Mormalized Temp. Grad. vs. Xdcer Position 
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Figure 19° Normalized temperature gradient as a function of transducer position 
required to produce 1.1 kW of total electric power. 





System Efficiency vs. Xdcr Position 
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Figure 20: System efficiency as a function of transduce: position. 





VI. CONCLUSIONS AND RECOMMENDATIONS 


It is not always obvious that the trends indicated by the figures in the proceeding 
chapter are those which we should expect. The lack of a working thermoacoustic 
generator for comparison does not allow for experimental validation of the resulis 
obtained. The results are of the correct order of magnitude, but the details of the 
relationships require additional investigation. It seems that a reasonable framework 
has been established for a model. 

One way in which the program may be tmproved is by inclusion of viscous and 
thermal effects in the resonator. If this were done then the transducer parameters could 
be adjusted so that its mechanical impedance would be equivalent to that of a rigid end. 
Comparison of the output of the program to existing experimental data would then be 
possible. This may serve to validate the output of the program to some degree. 

The next requirement ts that a thermoacoustic generator be built and tested. It 
need not be as large as the five inch engine at L.A.N.L. to measure general performance 
trends. If the program is modified as suggested above hen a much smaller device in 
which the viscous and thermal effects in the resonat | are substantial may be as easily 


modeled. A smaller device, working at a muci: !>.. pressure would not require such an 


expensive transducer. 
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